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Abstract 

Recent lattice QCD calculations, at physical pion masses and small lattice 
spacings that approach the continuum limit, have revealed that non-diagonal 
quark correlators above the critical temperature arc finite up to about 2 Tc- 
Since the transition from hadronic to free partonic degrees of freedom is merely 
an analytic cross-over, it is likely that, in the temperature regime between 
1-2 Tc, quark and gluon quasiparticlcs and prc-hadronic bound states can 
coexist. The correlator values, in comparison to PNJL model calculations 
beyond mean-field, indicate that at least part of the mixed phase resides in 
color-neutral bound states. A similar effect was postulated for the in-medium 
fragmentation process, i.e. for partons which do not thcrmalize with the 
system and thus constitute the non-equilibrium component of the particle 
emission spectrum from a deconfined plasma phase. Here, for the first time 
we investigate the likelihood of forming bound states also in the equilibrated, 
parton dominated phase above Tc which is described by lattice QCD. 
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1 Introduction 



The study of the QCD phase diagram and thermodynamics is receiving increasing 
attention in recent years. This field of physics is particularly appealing because the 
deconfined phase of QCD can be produced in the laboratory, in the ultrarelativistic 
heavy ion collision experiments at CERN SPS, BNL RHIC, CERN LHC and the 
future facilities at the GSI (FAIR) and in Dubna (NICA). On the other hand, 
lattice calculations on QCD thermodynamics are reaching unprecedented levels of 
accuracy, with simulations at the physical quark masses and several values of the 
lattice spacing approaching the continuum limit. In addition, the interpretation of 
lattice results in terms of phenomenological models is of fundamental importance in 
order to understand the microscopic structure of the QCD deconfined medium. The 
information that can be obtained from these complementary approaches will shed 
light on the features of QCD matter under extreme conditions, one of the major 
challenges of the physics of strong interactions. 

One fundamental question is the nature of the effective degrees of freedom in the 
temperature range 1-2 T^: the experimental results available so far show that the hot 
QCD matter produced in the laboratory exhibits robust collective flow phenomena, 
which are well and consistently described by relativistic hydrodynamics [HEllS]. The 
data are actually consistent with a near zero shear viscosity over entropy (tj/s) ratio, 
which signals the existence of a strongly interacting perfect fluid rather than a weakly 
coupled plasma state. In order to explain such a strong coupling between the degrees 
of freedom one either has to allow a strong enhancement in multi-parton interactions 
|1] or a modiflcation of partonic states, i.e. quasi-particles or color-neutral bound 
states. The quasi-particle model has been recently compared to the latest lattice 
QCD data [5] resulting in a large, temperature dependent, effective mass for quarks 
and gluons near Tc [HI El E] • However, the model cannot simultaneously reproduce 
bulk thermodynamics (pressure, energy density, interaction measure) and quark 
number susceptibilities [5l|9]. Alternatively, the interaction between partons might 
be enhanced due to the existence of a large number of binary bound states, mostly 
colored, in the QGP [10] . In addition, in-medium hadronization and the formation of 
color-neutral objects inside the partonic flreball, due to short formation times, have 
been postulated in Refs. [HI |T2] for the non-equilibrium component of the particle 
emission spectrum from a deconflned plasma phase. In this paper we address the 
possible presence of bound states also in the equihbrated component of the QGP 
phase. 

Correlations between non-diagonal quark flavors or between strangeness and 
baryon number have been proposed as a diagnostics to understand the nature of 
QCD matter immediately above the phase transition [13] . These observables can be 
calculated from flrst principles in lattice QCD, and several results have been pub- 
lished over the last few years [T^ IT5] . However, the rather large quark masses and 
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lattice spacings used in these simulations did not allow to draw a decisive conclusion 
so far. Very recently, new lattice results have become available for these observables, 
with simulations at the physical quark masses and finer lattice spacings approaching 
the continuum limit [161 IE]- These results, which are detailed in the next section, 
seem to favor a scenario in which bound states are present in the deconfined medium 
for a certain temperature range above T^. 

In this paper we are comparing the latest lattice QCD calculations of suscep- 
tibilities with the predictions of the Polyakov loop-extended Nambu Jona-Lasinio 
model [18]; the model correctly describes the chiral phase transition of QCD and 
incorporates aspects of the deconfinement phase transition through the coupling 
with the Polyakov loop. Our purpose is to understand the nature of the relevant 
degrees of freedom; we concentrate our analysis in the temperature region just above 
the phase transition. By going beyond the mean-field approximation and using the 
Monte Carlo method applied to the PNJL model |l9j, we incorporate fluctuations 
of the condensates (chiral, pion and kaon condensates) and of the Polyakov loop. 
In this way, we estimate quark number susceptibilities and baryon-strangeness cor- 
relations for a partonic quark-gluon plasma containing mesonic zero-modes and an 
additional degree of freedom (the Polyakov loop) which couples the different fla- 
vors. We attribute the difference with respect to lattice results to the presence of 
finite-momentum bound states in the QGP. 

2 Lattice QCD results 

Recent lattice calculations unambiguously show that the transition from the hadronic 
to the partonic system at zero baryo-chemical potential is an analytic cross-over [20] . 
In such cases the critical, or transition, temperature is determined by the inflection 
point in the temperature dependence of the relevant observables. The main quan- 
tities that are used to determine the transition temperature are the Polyakov loop, 
energy density and quark number susceptibilities for the deconfinement phase transi- 
tion, and the quark condensates for the chiral phase transition. The smooth behavior 
of all these QCD observables as functions of the temperature, which is evident in 
the most recent calculations that employ smaller lattice spacings and realistic quark 
masses, leads to interesting cross-over phenomena [21] • For example, in the left panel 
of Fig. [T] we show a comparison between the previously available lattice results for 
the light quark number susceptibilities [15] and the new Wuppertal-Budapest results 
obtained with physical quark masses and smaller lattice spacings [ITJ: it is evident 
that, for the most recent data, the transition is less steep. 

Both deconfinement, as expressed through the renormalized Polyakov loop or 
quark number susceptibilities, and chiral symmetry restoration, shown in the chiral 
condensates, experience therefore an extended transition region before reaching the 
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Figure 1: Left: Comparison between the lattice results for light quark number 
susceptibilities obtained with the stout [17], asqtad and p4 actions [12] (notice that 
on the horizontal axis the temperature is normalized by the different values 
obtained with the different actions). Right: comparison between the lattice results 
for light and strange quark number susceptibilities, obtained with the stout action 
at physical quark masses and Nt = 12 (from Ref. [I^). For the definition of A^^ and 



Nt see eq. (11). 



fully deconfined and chirally symmetric state. At any given temperature in this 
cross-over range these parameters could thus be interpreted as signalling a mixed 
phase of degrees of freedom where bound states or chirally broken states will co-exist 
with free quarks and gluons according to the relative values of the Polyakov loop or 
the chiral condensates. This transition region was called "semi-QGP" in Ref. [22], 
as opposed to the "full-QGP" at large temperatures, where the Polyakov loop is 
close to one and flat. 

Furthermore, the latest lattice results signal a flavor dependence of quark number 
susceptibilities even in the light quark sector. The rise of the strange quark suscep- 
tibility with temperature is slower and takes place at larger temperatures compared 
to the u case, as shown in the right panel of Fig. [Tj This feature was less pronounced 
in previous lattice results [15], since in that case nis/mu^d = 10, whereas the new 
results use a physical quark mass ratio {ms/mu^ = 28.15). 

This flavor difference between quark number susceptibilities in the light sector 
likely indicates that strange quarks experience deconfinement at slightly larger tem- 
peratures, compared to light quarks, thus implying a survival of strangeness-carrying 
hadrons in the QGP immediately above Tc. 

Another piece of evidence pointing in this direction is the behavior of the baryon- 
strangeness correlator Cbs or the non-diagonal quark number susceptibilities. Non- 
diagonal u — s susceptibilities (shown in Fig. [3]) exhibit a pronounced peak in the 
vicinity of the phase transition, and remain finite for relatively large temperatures 
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above T^. Similarly, the baryon-strangeness correlator (also shown in Fig. [3]) only 
reaches the predicted value for a purely partonic QGP (approaching the Stefan- 
Boltzman limit) near 2 T^. Although it was shown in Ref. [23j| that correlations 
between different flavors are nonzero in perturbative QCD at large temperatures 
due to the presence of flavor-mixing diagrams, the lattice data exhibit a strong 
enhancement of these correlations in the vicinity of T^, which survives up to relatively 
large temperatures above the transition [T7] and which cannot be accounted for by 
the perturbative QCD contribution alone. Taking into account this behavior, one 
could come to the conclusion that in the region 1-2 Tc the probability of forming 
color neutral bound states is quantifiable even in the case of a fully equihbrated 
system of quarks and gluons as simulated through lattice QCD. 

3 PNJL model 

The 2+1-flavor PNJL model is specified by the Euclidean action 

SEi^p, 0) = y dr J Sx [^^ dr^ + H{iP, ^\ 0)] + 5Se{(P, T) (1) 

with the fermionic Hamiltonian density 

•H = -27/^t(^. V + 74mo-0)z/' + V(V^,^^) , (2) 

where ip is the Nf = 3 quark field and rriQ = diag{mu,md,ms) is the quark mass 
matrix, with m„ = rrid 7^ m^. Quarks move in a background color gauge field 
= ^44 = iAf), where Aq = d^^gA'^t"- with the S'[/(3)c gauge fields and the 
generators = A'^/2. The matrix valued, constant field relates to the (traced) 
Polyakov loop $ as follows: 



Pexp ( i I drA^ 







^Tre^<^/^. (3) 



In a convenient gauge (the so-called Polyakov gauge), the matrix is given a diagonal 
representation 

= 03 A3 + 08 Ag , (4) 

which leaves only two independent variables, 03 and 08. The piece 5Se = —^U 
of the action ([T| controls the thermodynamics of the Polyakov loop. It will be 
specified later in terms of the effective potential, U{^,T), determined such that 
the thermodynamics of pure gauge lattice QCD is reproduced for T up to about 
twice the critical temperature for deconfinement. At much higher temperatures 

= 7o 7 and 74 = ijo in terms of the standard Dirac 7 matrices. 
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Parameters 



Physical quantities 



A [GeV] 


0.6023 




3.67 




24.72 


mou,d [MeV] 


5.5 


mo. [MeV] 


140.7 



U [MeV 


92.4 


\{^^)u,d\'^^ [MeV] 


241.9 


\{ip'^)s\^^^ [MeV] 


257.7 


[MeV] 


139.3 


tuk [MeV] 


497.7 



Table 1: Left: NJL model parameters. Right: physical quantities used to fix the 
parameters. 



where transverse gluons begin to dominate, the PNJL model is not supposed to be 
applicable. 

The interaction V in Eq. (j2]) is defined as follows: 



G 



V = -- ^ \{ipf^f) + {'^PfiJ5r^pfy 



f=u,d,s 



+ 



K 



det ^(1 + 75)V'i) + det ^(1 - 75)V^j) 



(5) 



Spontaneous chiral symmetry breaking is driven by the first term in Eq. ([s]), while 
the second term breaks the axial U{1)a symmetry explicitly. The NJL part of the 
model involves five parameters: the quark masses which we take equal for u- and 
(i-quarks and heavier for s quarks, the coupling strengths G and K and a three- 
momentum cutoff A. We take those from Ref. [21], they are listed in Table [l} The 
effective potential W($,T) which controls the dynamics of the Polyakov loop has 
the following form: 

T) = --a(T) + h{T) In [1 - 6 $*$ + 4 (<l>*^ + <l>=^) - 3 ($*$)^] (6) 
2 



with 

a(T) = ao + ai + 03 

The parameters are taken from the literature: 



h{T) 



T 



(7) 



Co = 3.51 , ai = -2.47 



a2 



15.22 , 6, = -1.75 



The critical temperature Tq for deconfinement in the pure gauge sector is fixed at 
270 MeV in agreement with lattice results. 

In the model, quarks acquire a constituent mass through their interaction with 
the chiral condensate. Due to the fiavor-mixing term in the Lagrangian, the mass of 
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a given flavor also gets contributions from the chiral condensates of the other quark 
flavors: 



rrii = rriQi 



K 



Wj){<^k) = rrioi - 2G{i)ii)i) + K{ipj^pj)l(il)ki)k) i ^ j ^ k. 



In the PNJL model the partition function in momentum space is written as 



/ 



V(f)T>(Ta VlTa'DSa 'DPa CXp 



V / 1 



T\2 ^1 (27r)3 



lndet[S' ^{iun,p,fiq)] 



where Un = (2n + 1)ttT are the Matsubara frequencies, cTq, and tTq, are 18 bosonic 
fields corresponding to the possible scalar and pseudoscalar condensates while 
and Pa are additional 18 auxiliary fields necessary in order to deal with the six- 
fermion interaction induced by the 't Hooft term. For symmetry reason, clearly 
(tTq,) = at the mean field level, but we will let these fields fluctuate in the Monte 
Carlo approach, too. This will allow to take the contribution of mesonic zero- 
modes into account. ^^(0, T) is the Polyakov loop potential given above and the 
S(aa, Ha, Sa,Pa) IS the bosonic action 



obos 



G 

(JaSa + T^aPa + \SaSa + PaPa\ + 



K 

+ ^^a/37 [SaSpS^ — diSaPpP^y] 



(9) 



where the constants Aap-y are expressed in terms of the Gell-Mann matrices accord- 
ing to 



A/37 '■= -^fijkemni {K)im {^p) jn i^j) kl «, /3, 7 G {0, . . . , 8}. 

In order to perform the integration in Sa and Pa, the stationary phase approxi- 
mation is used, choosing the fields Sa and Pa so as to minimize the integrand in the 
bosonized partition function. The necessary condition imposed on the fields is 



3K 

Oa + GSa H —Aal3^ ['S'/^S'^ — -P3-P7] = , 

3K 

7Ta + GPa —Aap-yS/sPj = . 



(10) 



The presence of a volume factor V in the exponent of Eq. [8] makes it possible 
to compute the full partition function using Monte-Carlo techniques. In this way 
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we consider not only the saddle-point contributions, but also configurations that 
correspond to fluctuations around the minima of the action. 

The size of the volume is now specified according to the conventions adopted 
in lattice calculations. For a fixed extension of the lattice in the Euclidean time 
direction, the temperature is set by the lattice spacing a, and the volume size is 
related to the temperature: 



1 



a — 



(11) 



where Nt is the number of lattice sites in the Euclidean time direction, and Ns is 
the number of lattice sites in the space direction. It follows that 



V = k/T\ 



(12) 



In particular here we used the typical values of Ng/Nt used in Lattice simulations 
which correspond to A; = 27 and k — 64. 

Let us see now the definition of quark number susceptibilities. The thermody- 
namic potential can be expanded in a Taylor series in ^iq/T around zero chemical 
potential. 



with 



n(T,rt = ^i5'"2=EC(r)('^)*(Y)' 



1 



i,j=0 



d{fim/Tyd{iin/Ty 



=Mn=0 



(13) 



(14) 



where only even terms survive due to CP symmetry. The coefficients c^"'{T) are 
evaluated at Hq — 0. The baryon-strangeness correlator is defined in the following 
way: 



C 



BS 



1 + 



ds 



(15) 



Looking at the definition of the partition function we immediately see that the 
susceptibilities c^"' involve derivatives of lndet[S~^{iu!n,p, fJ^q)], which is the only 
term in Z which explicitly depends on the chemical potentials. 
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For the first derivative of the log of the partition function we have: 
d]nZ{T,iiu,IJ'd,IJ's) ^ 

In / Da^r'Tr^r'^exp — Indet /i^, /i^) - >Sg[^ 



1 V 
V /dlndet S~'^{fiu, fJ^d, fJ^^ 



(16) 



Going on in the same way for the second derivative, we obtain the coefficients cf^'^ 
(quark number susceptibilities) 

1 92 72/^/^2 

cfi'^' = — — In^: = — — — ( — IndetS" ^{^u, jJ^d, jJ^s) 



VT dpig, dpi,, VTHT\ dpi,, dpi, 



(^) ((9^^^'^^^'^ \l^u, IJ'd, IJ's)^ ^^^Indet^" \l^u, IJ'd, IJ's)^^ 

lndetS~^{liu,IJ'd,IJ's)^ ^^^lndet5'"^(/x„,/Xd,/^s)^ j • (17) 



VV I d 



Therefore, the form of the fermionic determinant and how we introduce the cutoff 
are crucial for our calculation. 

Our choice will be the following: we consider the effect of the condensates only 
for momenta smaller than the cutoff A, while in the high p region only free quarks 
are included in the calculation, namely {M{aa, tTo) = 'S'~^((Ta, tTq,)): 

V- P d^P 1 1 X V- /"^ d^P 1 1 

/ 77r-T^lndetM(cr„,7ra,a;n) = >^ / 77r^lndetM(cT„,7ra,a;n) 
^ Jo (27rj^ [27ry 

+ Z2j^ ^^lndetM(a„ = 0,7ra = 0,a;„). (18) 

For the free term of the decomposition we know the eigenvalues of the fermionic 
matrix analytically and therefore the sum over the Matsubara frequencies can be 
performed exactly. In this way the fermionic term can be rewritten as: 

/ TTT^ In det M((Ja, TTa, U;n) = / —--Indet M{aa,'Ka,UJn) 
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Figure 2: Renormalized chiral condensate A; g: comparison between the lattice re- 
sults (light blue band) [21], the mean-field PNJL model result (full, black line), and 
the Monte-Carlo PNJL model result (blue, dash-dotted line). 

where the sum over the index j means the sum over the different eigenvalues of the 
fermionic matrix. 

The importance of including fiuctuations in the model is shown in Fig. |2j where 
we compare the chiral condensate from lattice QCD (from Ref. |2l]), the PNJL 
result at the mean-field level, and the PNJL result in which fiuctuations of all fields 
are included. As it is evident, the inclusion of fiuctuations makes the curve much 
smoother and it brings it to a good agreement with the lattice data (a similar effect 
was observed in Ref. [25]). 

4 Results and conclusions 

In order to draw conclusions on the presence of bound states in the QGP, the most 
relevant comparison between PNJL and lattice results is for the non-diagonal quark 
correlators. In particular, we will focus here on the u — s correlator and on the 
baryon-strangeness correlator, for which new lattice results have been recently pub- 
lished [17j. At the mean field level, the PNJL model has no correlations between 
the different quark flavors, therefore the corresponding u — s correlator stays flat 
and equal to zero over the full temperature range. Similarly, the baryon-strangeness 
correlator takes everywhere the value corresponding to the non-interacting QGP, 
namely one [13] . In order to properly estimate all possible contributions to this ob- 
servable from colored degrees of freedom and mesonic zero modes, we need to go be- 
yond mean fleld and take fluctuations of all flelds (Polyakov loop, chiral condensates, 
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Figure 3: Left: Comparison between the lattice results for the u — s correlator as a 
function of T/Tc [17_\, and the PNJL model results. The mean field PNJL result is 
zero for all temperatures, as expected (dashed curve). The blue curve corresponds 
to the PNJL model result when only the Polyakov loop fluctuations are taken into 
account. The red curve is the full PNJL model prediction, with fluctuations of all 
fields taken into account. Notice that the red curve will fall on the blue curve in the 
infinite volume limit. Right: Baryon-strangeness correlator: comparison between 
the continuum extrapolated lattice results from Ref. [17] (red band), the PNJL 
model result at the mean field level (black, dashed line) and the PNJL model result 
with inclusion of Polyakov loop fluctuations (blue, solid line). 



pion and kaon condensates) into account. By plotting the results corresponding to 
fluctuations of the Polyakov loop only, we will be able to determine the relative 
strength of the quark correlator due to colored states above Tc. The difference be- 
tween such a PNJL calculation and lattice QCD can therefore give us an estimate on 
the relative abundance of bound states in the medium above T^. The fluctuations 
of the condensates give us an estimate of the contribution due to the zero-mode 
mesonic states. We will find that these fluctuations vanish in the infinite volume 
limit, leaving the Polyakov loop fluctuations as the only non-zero contribution to 
this observable. The results of our calculation are shown in Fig. [3| since the PNJL 
model describes the deconfined state, we only show the corresponding curves for 
T > Tc. In the left panel, the lattice results for the u — s correlator as a function of 
T/Tc are compared to the PNJL model results at the mean field level (dashed line), 
when Polyakov loop fluctuations are taken into account (blue line) and when fluctu- 
ations of all fields are taken into account (red line). The same comparison is made 
in the right panel for the other observable we are considering, namely the baryon- 
strangeness correlator. It is evident from the left panel, that even a PNJL model 
taking into account all possible fluctuations cannot fully account for the strong dip 
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and the slow rise of the non-diagonal quark correlator determined by lattice QCD. 
Notice that the red curve has been obtained in a Monte Carlo simulation for which 
the ratio Ng/Nt (see Eq. (11)) is the same as the one used in the lattice simulations. 
However, we find that in the thermodynamic limit of Ns — ?■ oo, the blue curve, 
which corresponds to the colored QGP contribution to this observable, remains the 
same, while the red curve falls on top of the blue curve, i.e. fluctuations of the 
condensates (zero-mode mesonic contributions) vanish in the infinite volume limit. 
It is evident that the lattice data can be reproduced only for large temperatures, 
thus allowing for a considerable contribution from finite-momentum bound states of 
both baryonic and mesonic nature. 

In the right panel we show the baryon-strangeness correlator. A comparison 
between lattice data and PNJL model results for this observable again suggests the 
presence of bound states in the QGP for temperatures up to 1.6-1.7 Tc. 

Furthermore, the peculiar shape of the non-diagonal correlator near Tc (sharp 
dip and subsequent slow rise towards zero) can be interpreted when comparing it 
to the separate contributions of mesonic and baryonic states in a hadron resonance 
gas (HRG) calculation [26]. Mesonic states in the HRG model exhibit a negative 
correlation, whereas baryonic states yield a positive value (see Fig. |4]). The dip and 
slow rise thus is likely caused by enhanced baryonic state formation (or survival) 
at higher temperatures. Still, the lattice data never exceed zero before full decon- 
finement is reached, which means that baryonic states never dominate the hadron 
formation. This is confirmed by the magnitude of the difference between PNJL and 
lattice QCD in Cbs- 

Taking into account the flavor dependence of the light quark susceptibilities as 
shown in the right panel of Fig. [T| one can deduce a scenario where strange quark 
bound states are formed (or survive) at higher T in the deconfined medium than 
light quark bound states. Earlier lattice calculations only exhibited this effect in the 
comparison between light and heavy quarks, but the recent improvement in lattice 
accuracy indicates effects already at the strange quark level, which leads to specific 
experimentally verifiable predictions, such as an enhanced survival probability of 
strange over non-strange resonances near, but above, Tc. 

In summary, a comparison of PNJL and lattice QCD calculations yields ample 
evidence, that a phase of mixed degrees of freedom exists for a particular temperature 
range above the QCD transition temperature. The PNJL model calculations and the 
mapping of the flavor and baryon number dependence in lattice QCD calculations, 
show that it is likely that in this equilibrated phase the relative abundance of strange 
hadronic and baryonic states is larger at high T rather than immediately above Tc, 
but mesonic bound states still dominate over the entire temperature range. In 
other words, strange hadrons and non-strange baryons form earlier in the mixed 
phase than their non-strange and mesonic partners. These dependencies should 
be experimentally verifiable through particle identified measurements of light- and 
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Figure 4: u — s correlator: the different curves show the contributions of baryons 
(blue), mesons (red) and the total (black) from the Hadron Resonance Gas model. 

strange-quark hadrons and resonances near the phase boundary. 
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